Ballistic magnon transport and phonon scattering in the antiferromagnet Nd 2 Cu04 
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The thermal conductivity of the antiferromagnet NcbCuCU was measured down to 50 mK. Using 
the spin-flop transition to switch on and off the acoustic Nd magnons, we can reliably separate 
the magnon and phonon contributions to heat transport. We find that magnons travel ballistically 
below 0.5 K, with a thermal conductivity growing as T 3 , from which we extract their velocity. We 
show that the rate of scattering of acoustic magnons by phonons grows as T 3 , and the scattering of 
phonons by magnons peaks at twice the average Nd magnon frequency. 

PACS numbers: 72.15.Eb, 74.72.-h, 75.30.Ds, 75.50.Ee 



The collective excitations of an antiferromagnet are 
magnons. Like electrons and phonons, magnons can 
transport heat. Several groups have recently used heat 
transport toprobe spin excitations in quantum magnets 
0, Sll II U @ The results tend to be difficult to in- 
terpret because magnons interact strongly with phonons. 
However, as T — > 0, the two become decoupled and their 
respective mean free paths are expected to reach the size 
of the sample. This ballistic regime is easily reached 
for phonons, but it has rarely been demonstrated for 
magnons, because acoustic magnons are almost always 
gapped. If they were not gapped, their thermal conduc- 
tivity would be hard to separate from that of acoustic 
phonons, as both are expected to exhibit a T 3 tempera- 
ture dependence in the ballistic regime. In this context, 
the antiferromagnetic insulator Nd2Cu04 offers a rare 
opportunity. Jin et al. have recently shown that in this 
material magnon heat transport can be switched on and 
off with a magnetic field via a spin-flop transition • 

In this Letter, we use this "switch" to investigate the 
physics of magnons. First, by cooling down to 50 mK, 
we show that the ballistic regime for magnon transport, 
characterized by a clear T 3 conductivity, sets in below 0.5 
K. The magnitude gives us the magnon velocity, whose 
value agrees with a recent calculation |8|. Secondly, we 
show how above 0.5 K magnons are scattered by phonons 
at a rate proportional to T 3 . Finally, we show that the 
scattering of phonons by magnons is strongly peaked at 
a temperature where the optimum phonon energy equals 
twice the average magnon energy, in line with a model of 
dominant one phonon- two magnon processes [9J. 

Nd2Cu04 is the insulating parent compound of the 
electron-doped high-T c superconductor. The spins on the 
Cu 2+ ions order antiferromagnetically below Tjy « 275 
K, with spins lying in the Cu02 planes. As the temper- 
ature is lowered, the ordered spin structure undergoes a 
series of re-arrangements, to reach a non-collinear state 
at low temperature (T < 30 K), with the spin direc- 
tion chan ging by 90 degrees in alternate planes along the 
c-axis [sl lint Flit ll2T ]. Due to the magnetic exchange 



interaction between Cu 2+ and Nd 3+ , moments on the 
Nd site also order in the same structure. The magnon 
spectrum of Nd2CuC>4, studied by neutron scattering 
HH Qja LL-il UM > is characterized by two sets of branches, 
associated respectively with spins on Cu and Nd sites. 
Cu magnons have energies above 5 meV 0, and 
four optical Nd magnon branches lie in the range 0.2 to 
0.8 meV 0IHJ. A recent model predicts that the Nd 
spin sublattice should also support an acoustic magnon 
mode, with a slope of 10 meV A and a very small gap 
(5 //eV) |8j. An in-plane magnetic field of 4.5 T causes a 
spin-flop transition from non-collinear to collinear state. 
This change of spin structure modifies the Nd magnon 
spectrum, as evidenced by the additional heat conduc- 
tion which appears at the transition . 

Single crystals of Nd2CuC>4 were grown by a standard 
flux method. Our sample was annealed in helium for 10 
h at 900 °C, and cut to a rectangular shape of dimensions 
1.50 x 0.77 x 0.060 mm 3 (length x width x thickness in 
the c direction). Contacts were made with silver epoxy, 
diffused at 500 °C for 1 hour. The thermal conductivity 
k(T) was measured using a standard steady-state tech- 
nique |l6j. in a dilution refrigerator (50 mK to 2 K, with 
R11O2 thermometers), and in a 4 He cryostat (above 2 K, 
with Cernox thermometers). The magnetic field H was 
applied perpendicular to the heat current, both in-plane 
and along the c-axis. Note that in a magnetic insula- 
tor, the thermal conductivity is the sum of phonon and 
magnon conductivities: k(T) = k p (T) + K m (T). 

Ballistic regime. - Fig. 1 shows the temperature 
dependence of k, in H = 0. 10 T 1 c. and 10 T // c [l7|. 
The in-plane magnetic field causes a huge enhancement 
of discovered by Jin et al. above 2 K 0]. The 

advantage of going lower in temperature by a factor 40 is 
that one can then unambiguously decouple k p and n m . In 
Fig. 2 we plot the zero-field data (upper panel), as «(0) 
vs T 26 , and the difference between in- field data and zero- 
field data (lower panel), as An — k(10 Tic)- k(0) vs 
T 3 . In both cases, the data follows a straight line below 
about 0.5 K, which defines the ballistic regime for both 
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FIG. 1: Temperature dependence of the in-plane thermal 
conductivity k of Nd2Cu04, in zero magnetic field (circles; H 
= 0) and in a field of 10 T applied either in the basal plane 
(triangles; H _L c) or along the c-axis (squares; H/ /c). Inset: 
Zoom below 1 K. 



types of carriers, wherein the mean free path is limited 
by sample boundaries. 

a) Phonons (T < 0.5 K). As T — ► 0, the phonon mean 
free path in insulators becomes limited only by the phys- 
ical dimensions of the sample. In this ballistic regime, 
the phonon conductivity is given by [3: 



k p = —-K 2 k B (k B T/Kf <v p 2 > l , 



(1) 



where <v~ 2 > is the inverse square of the sound velocity 
averaged over the three acoustic branches in all q direc- 
tions, and Iq is the temperature-independent mean free 
path, given by the cross-sectional area A of the sample: 
Iq = 2yJ A/n. For our sample, Iq = 0.24 mm. 

For samples with smooth surfaces, a fraction of 
phonons undergo specular reflection and the mean free 
path exceeds Iq, estimated for diffuse scattering. Because 
specular reflection depends on the ratio of phonon wave- 
length to surface roughness, the mean free path depends 
on temperature and k p no longer grows strictly as T 3 
181. The general result is that n p oc T a with a < 3 
19| . Typical values of a are 2.6-2.8 (see [l6| and refer- 
ences therein) . For instance, n p oc T 2 7 in the insulator 
La2CuC>4 j2fj (where all magnons are gapped). In Fig. 
2, the zero- field conductivity of our Nd2CuC>4 crystal can 
be fitted to k p = bT a below 0.4 K, with a = 2.6. This 
is a first indication that k(0) is mostly due to phonons, 
with a ballistic regime below ~ 0.5 K. 

The second argument comes from the magnitude of b. 
Taking Iq — 0.24 mm as a lower bound on the mean free 
path and applying Eq. (1) at T = 0.4 K, where k(0) 
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2: Top panel: Zero-field thermal conductivity k(0) vs 
The heat at H — is carried entirely by phonons: n p 
~ s(0) (see text). The line is a linear fit below 0.4 K. Bottom 
panel: The additional thermal conductivity induced by a 10 
T in-plane field, Aac = k(10 Tic)- «:(0), plotted as Ak vs 
T 3 . The line is a linear fit below 0.5 K. The additional heat 
is carried entirely by acoustic magnons: K m = Ak (see text). 



= 0.61 mW / K cm, we obtain: v p = (< v p 2 >) -1 / 2 
= 3.2 x 10 3 m/s. A mean free path larger by a factor 
of 2 would give v p — 4.5 x 10 3 m/s. These values are 
in good agreement with the measured sound velocities of 
Nd 2 Cu0 4 : v L = 6.05 x 10 3 m/s and v T = 2.46 x 10 3 m/s 
for longitudinal and transverse modes along [100] fij. So 
it appears that phonons account for the full k(0). In a 
third argument presented later, we make the case that 
magnon transport is negligible in the non-collinear state 
at any T, so that k(0) ~ n p at all T. 

b) Magnons (T < 0.5 K). Magnon transport is re- 
vealed by looking at the difference between the two curves 
in the inset of Fig. 1, Ak = k(10 Tic)- k(0), plot- 
ted in the bottom panel of Fig. 2. We get the strik- 
ing result that Ak — dT 3 , below 0.5 K. Such a field- 
induced additional T 3 contribution to k at low T can be 
attributed unambiguously to acoustic magnons. To our 
knowledge, this is the first observation of the characteris- 
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tic T 3 dependence of magnon thermal conductivity |22fl . 
As bosons with a linear energy dispersion at low q, acous- 
tic magnons should contribute to heat transport in the 
same way as acoustic phonons, provided there is no gap 
in their spectrum. Hence for gapless acoustic magnons, 
Eq. 1 should apply in the ballistic regime, with <v~ 2 > 
replacing <v~ 2 > [23]. Note, however, that the prefac- 
tor will depend on the number of magnon branches. In 
the simplest case, one expects two degenerate magnon 
branches, with the degeneracy lifted by a magnetic field. 
As our data is in 10 T, we assume that only one branch 
lies at zero energy. Hence: 

% = rx ^-n 2 k B {k B T/h) 3 <v m 2 > l . (2) 
o 15 

Assuming that k(0) = k p and that k p is independent of 
magnetic field in the ballistic regime, we have n rn = k 
- k p = Aft, so that K m = dT 3 , with d = 15 mW K~ 4 
cm . Using Eq. 2 with — 0.24 mm, we get v m = 
(<w m 2 >)- 1 / 2 = 1.5 x 10 3 m/s. 

With the important caveat that their calculation was 
done for the non-collinear state (H = 0), our findings 
agree with the theoretical calculations of Sachidanandam 
et al. @: 1) the acoustic mode they predicted is con- 
firmed, 2) the predicted magnitude of the gap (5 ^eV) is 
consistent with our lowest temperature (50 mK), and 3) 
the mode dispersion (or magnon velocity) is of the right 
mag nitude (10 meV A = 1500 m/s). 

An intriguing difference between phonons and 
magnons in the ballistic regime is that only the former 
appear to undergo specular reflection. Indeed, the per- 
fect T 3 dependence of n m implies that the fraction of 
magnons that are specularly reflected must be negligible, 
while it is not for phonons. We speculate that the strong 
anisotropy expected of the magnon velocity in Nd2CuC>4, 
compared to the roughly isotropic phonon velocity, might 
account for this difference. It is also likely that the spin 
structure is disordered at the surface so that the surface 
is rougher magnetically than structurally. 

Inelastic regime. — Beyond the ballistic regime, the 
temperature interval up to 10 K or so is dominated by 
phonon- magnon scattering: a) acoustic phonons are scat- 
tered strongly by magnons, and b) acoustic magnons are 
scattered strongly by phonons. 

a) Phonons (T > 0.5 K). A 10 T field along the c-axis 
leaves the system in its non-collinear state and as seen 
in Fig. 1, it has negligible impact below 1 K. Above 
1 K, however, it clearly suppresses k, relative to k(0). 
The 10 T // c curve is striking: the rapid growth in 
phonon conductivity generally observed in insulators is 
almost entirely stopped between 1 and 3 K. An extremely 
effective inelastic scattering mechanism switches on near 
1 K and off above 3 K. This range of energies is precisely 
that of Nd magnons in NxbCuO,}. 

The theory of phonons scattered by magnons was given 
a thorough treatment by Dixon |9j , who argues that the 
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FIG. 3: Top panel: Normalized thermal magneto-resistance 
vs temperature, for H//c. This is the field-dependent thermal 
resistivity of phonons. The arrow indicates the temperature 
at which magnon scattering is maximum: T* = 3 K. T* is 
directly related to the magnon energy < E magrl on > via the 
peak phonon frequency, uip honon (see text). Bottom panel: 
Difference in thermal conductivity caused by rotating a 10 T 
field from in-plane to c-axis, An' — k(10 Tic) - k(10 T // 
c), vs temperature. This is the conductivity of magnons. The 
dotted line is the estimated magnon specific heat (see text) 
and the solid line is a fit as indicated. 



one phonon- two magnon process is dominant. Dixon gets 
two main results: 1) the scattering peaks at a tempera- 
ture such that the optimum phonon frequency (~ Ak B T) 
is equal to the average energy of two magnons, and 2) a 
magnetic field enhances the effectiveness of magnon scat- 
tering, especially when the field is perpendicular to the 
spin direction. We attribute the suppression of k induced 
by applying H//c to this enhanced effectiveness. 

Phonons are scattered by magnons and several other 
processes (boundaries, defects, other phonons). However, 
because only magnon scattering depends on magnetic 
field, we can bring out the role of magnon scattering by 
focusing on the change in thermal resistivity, w = l/«, 
i.e. the thermal magnetoresistance Aw = w(10 T // c) - 
w(0). In the top panel of Fig. 3, the normalized magne- 
toresistance, Aw/w(0), is plotted vs T up to 20 K. Quite 
evidently, phonon transport emerges as a spectroscopy 
of magnons. (The same is true of electron heat trans- 
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port, as was recently shown for CeRhIn 5 H3-) Indeed, 
the maximum scattering of phonons will occur when the 
peak phonon energy matches twice the average energy 
of optical magnons, < E magnon > 0- At T = T*, 
heat transport is dominated by thermal phonons with 
frequencies around u* honon — AksT* /h. Therefore, in 
Nd2CuC>4, where < E magnon > = 0.5 meV is the average 
of the measured optical magnon energies (which range 
from 0.2 to 0.8 meV 0E3), w(T) is expected to peak 
when hw* honon = 2 < E magnon >~ Ak B T*, i.e. k B T* 
0.2") meV, or T* ~ 3 K, precisely as observed. 

b) Magnons (T > 0.5 K). We assume that magnons in 
10 T // c scatter phonons just as well as in 10 T ± c, 
so that « p (10 T // c) ~ Kp(10 T _L c), and « m (10 T 
_L c) ~ «(10 Tic)- k(10 T // c) = Ak' (which is 
equal to Ak below 1 K). In the bottom panel of Fig. 3, 
we plot Ak' vs T. The magnon conductivity n m peaks 
at 3 K, because the magnon heat capacity C m (T) satu- 
rates when T exceeds the maximum energy of the acous- 
tic magnon branch, while the inelastic scattering due to 
phonons keeps increasing. That maximum energy is un- 
known, but we expect it to lie roughly in the range 0.2- 
0.8 meV. We estimate C m (T) using an Einstein model: 
C m = Ce = z 2 e z (e z -1)~ 2 where z = hujE/ksT, with 
Huje = 0.5 meV. This is clearly inaccurate at very low T 
where C m oc T 3 , but should be reasonable for temper- 
atures comparable to and above the maximum energy, 
given roughly by Hue- As seen in Fig. 3, C m is con- 
stant above 5 K or so and the temperature dependence 
of K m comes entirely from the scattering rate, which we 
model simply as the sum of boundary and phonon scat- 
tering: r = To + T p . As seen in Fig. 3, a good fit to 
K m (T) is obtained with T p oc T 3 (and Fq constant, as it 
should) . This shows that acoustic phonons scatter acous- 
tic magnons (as they do electrons) in proportion to their 
number, which increases as T 3 . The magnon mean free 
path in Nd2Cu04 decreases by a factor of 15 between 5 
and 20 K, and should continue to decrease up to roughly 
the Debye temperature. This suggests that the magnon 
mean free path is unlikely to ever be independent of tem- 
perature in the range 20 to 200 K, as was assumed in the 
analysis of data on La2Cu04 [jfl. 

In summary, by cooling to 50 mK, the characteristic 
T 3 heat transport expected of ballistic magnons has for 
the first time been observed, revealing the existence of 
an acoustic magnon branch (gapless or with a very small 
gap, with upper bound ~ 0.1 K) in the collinear phase of 
Nd2CuC>4 (H > 7 T _!_ c). Paradoxically, the properties 
of this mode agree well with recent calculations applied 
to the non- collinear phase (H = 0) 8], where our data (in 
H = or H//c) rules it out. A calculation specifically for 
the collinear phase in field is desirable, as well as a search 
for the acoustic branch via low-energy neutron scatter- 
ing. Above the ballistic regime, phonons scatter magnons 



strongly, initially at a rate T p oc T 3 , so that the magnon 
mean free path decreases rapidly with temperature, and 
should continue to do so up to roughly the Debye temper- 
ature. Phonons are themselves strongly scattered by the 
magnons, in the temperature range where one phonon- 
two magnon processes are most effective, namely in the 
vicinity of 4fceT = 2 < E magnon >. In general, one 
can expect phonons in all cuprates (at any doping) to be 
scattered significantly by magnetic excitations (typically 
from Cu moments). 
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